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ABSTRACT 

The temperature dependences of the lifetime of polycubanes C4+4 n H 8 with n = 2 - 
5 up to their decomposition have been directly calculated using the molecular dynam- 
ics method. It has been shown that the activation energy of decomposition of these 
metastable clusters, in which the C-C bonds form an angle of 90° that is not characteris- 
tic of carbon systems, rapidly decreases with an increase in n due to the lowering of the 
energy barrier that prevents the decomposition of the clusters. This has cast some doubt 
on the recently made suggestion that there exist nanotubes {n » 1) with a square cross 
section. Nonetheless, the stability of bicubane {n = 2) and tricubane {n = 3) has proved 
to be sufficient for their existence at the liquid-nitrogen temperature. 



1. INTRODUCTION 

In the best known carbon compounds (diamond, graphene, carbyne), the angles between 
covalent bonds are equal to 109.5°, 120°, and 180°, which correspond to the sp 3 , sp 2 , and 
sp hybridizations of atomic orbitals, respectively (in fullerenes and carbon nanotubes, 
these angles differ only insignificantly from 120°). Consequently, increased interest has 
been expressed by researchers both in the C 8 H 8 cubane (Fig. 1), which was synthesized 
for the first time in 1964 [1] and where the angles between all adjacent C-C-C bonds are 
equal to 90°, and in the supercubane (Fig. 2) consisting of C 8 "cubes" covalently bound 
to each other in the directions of their main diagonals (note that, although there are 
publications on the synthesis of the supercubane, the fact of its existence, as far as we 
know, has still not been proved [2, 3]). Since so strong "bending" of the C-C-C bonds is 
energetically unfavorable, the cubane and the supercubane have a substantially strained 
structure. The cubane structure is stabilized by hydrogen atoms (Fig. 1), whereas the 
stability of the supercubane is provided by the interaction of neighboring elements C 8 
(Fig. 2): in both cases, there is a saturation of "dangling" bonds of the carbon atoms at 
the vertices of the cube (an isolated cubic cluster C 8 is unstable [4]). 

The question now arises as to the possibility of the existence of pure carbon structures 
that do not contain hydrogen atoms (except, perhaps, for the passivation along the edges) 
and in which all C-C bonds are either perpendicular or parallel to each other (in the 
supercubane, this is true only for the C-C bonds within each individual cubic element). 
In [5], the polycubanes C 4+4n H 8 with n = 2 - 6 (Fig. 3) were theoretically investigated 
using the density functional theory method [For polycubanes, the authors of [5] used the 
structural formula C 8+8n H 8 (n > 1). We believe that the formula C 4+ 4„H 8 (n > 2) is 
more convenient, because, in this case, n is nothing less than the number of C 8 cubes in 



the carbon skeleton, which have (at n » 1) common faces, so that n — 1 corresponds to 
cubane, and n — 2, 3, 4 , to bicubane, tricubane, tetracubane, etc.]. It was shown that all 
their vibration frequencies are real; i.e., these hydrocarbon clusters represent metastable 
(corresponding to local minima of the potential energy) atomic configurations. Later [6], 
it was noted that, since an increase in n leads to a decrease in the relative number of 
hydrogen atoms in the polycubane, at n >> 1 we are actually dealing with a thinnest 
carbon nanotube that has a square cross section (hence, here, the term "nanobar" is more 
appropriate). According to [6], this nanotube can be produced under the deformation of 
a conventional (2, 2) nanotube. 

It should be noted, however, that, although the real values of the vibrational frequen- 
cies of the system provide a necessary condition for its metastability, the quantitative 
characteristic of the degree of stability of the metastable state is the activation energy of 
its decay E a , which is involved in the Arrhenius formula for the temperature dependence 
of the lifetime r of the system in this state; that is, 



where A is the frequency factor with dimensionality s^ 1 , T is the temperature, and k B 
is the Boltzmann constant. As a rule (even though not always), the value of E a is close 
to the height U of the minimum energy barrier separating this metastable state from the 
configurations with a lower energy. Previously [7] , using the molecular dynamics method 
("computer experiment"), we directly calculated the dependence of the lifetime r on the 
temperature T for the C 8 H 8 cubane and demonstrated that the activation energy of its 
decomposition is E a = 1.9 ± 0.1 eV. Such a large value of E a indicates a high thermal 
stability of the cubane and explains why the cubane molecules not only retain their 
structure at temperatures significantly higher than room temperature [8], but also can 
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form a molecular crystal, the so-called solid cubane s-C 8 H 8 with a melting temperature of 
approximately 400 K [9] (during melting of this compound, only the weak van der Waals 
bonds between the C 8 H 8 clusters are broken, but the clusters themselves retaine their 
structure). 

The main objective of this work is to perform the numerical simulation of the dynamics 
of polycubanes C4 + 4„H 8 with n > 2 over a wide temperature range, as well as to determine 
the dependences of the activation energies of their decomposition E a and frequency factors 
A on the number n. We have restricted ourselves to the first four (following the cubane) 
representatives of the family: bicubane, tricubane, tetracubane, and pentacubane (n = 
2, 3, 4, and 5, respectively). This is enough to make sure that the thermal stability of 
the polycubanes drastically decreases with increasing n. We have also investigated the 
decomposition of all the aforementioned hydrocarbon clusters, determined the heights U 
of the energy barriers preventing their decomposition, and demonstrated that, within the 
limits of error, they coincide with their corresponding values of E a . This paper is organized 
as follows. Section 2 contains the description of the computational methods. Section 3 
presents the results (including the data obtained from the molecular dynamics calculation 
and static simulation) and their discussion. Section 4 provides the brief conclusion and 
inferences. 



2. COMPUTATIONAL METHODS 

The energies of arbitrary atomic configurations were calculated within the framework 
of the nonorthogonal tight-binding model [10], which represents a reasonable compromise 
between more rigorous ab initio approaches and extremely simplified classical potentials of 
the interatomic interactions. This model adequately describes the structural and energy 



characteristics of various hydrocarbon systems, including both small clusters (cubane 
[7], methylcubane [11], etc.), and macroscopic systems (graphane [12], graphone [13], 
hydrogen interstitials in diamonds [10], etc.). Yielding to the accuracy of the ab initio 
methods, this model requires a considerably less expenditure of computer resources and, 
therefore, makes it possible to investigate the evolution of the system consisting of 10 - 
100 atoms for a time (from 1 ns to 1 /is) sufficient to recruit the necessary statistics [14]. 
Previously, we, in particular, successfully used it to simulate the spontaneous regeneration 
of the graphene/graphane interface during thermal disordering [15], to determine the 
dependence of the dielectric gap of graphane nanoribbons from their width [16], etc. 

For the simulation of the dynamics of polycubanes, at the initial instant of time, random 
velocities and displacements were imparted to each of the atoms in such a manner that 
the momentum and the angular momentum of the cluster as a whole were equal to zero. 
Then, the forces acting on the atoms were calculated using the tight-binding model and 
the Hellmann-Feynman theorem [10]. The classical Newton equations of motion were nu- 
merically integrated using the velocity Verlet method with the time step t = 2.72 x 10~ 16 
s [14]. In the course of the simulation, the total energy of the cluster (the sum of the po- 
tential and kinetic energies) remained unchanged, which corresponds to a microcanonical 
ensemble (the system is not in a thermal equilibrium with the environment [7, 11-15]). 
In this case, the "dynamic temperature" T is a measure of the energy of relative motion 
of the atoms and is calculated from the formula [17, 18] (E^ n ) = ^k B T(3n — 6), where 
(i?kin) is the time-averaged kinetic energy of the cluster and n is the number of atoms in 
the cluster (here, it is taken into account that the cluster as a whole neither moves nor 
rotates; therefore, the number of degrees of freedom is reduced by six). 

In order to determine the heights U of the energy barriers preventing the decomposition 



of the polycubanes, we investigated the hypersurface of the potential energy E pot for each 
value of n as a function of the atomic coordinates in the vicinity of the initial metastable 
configuration and used the method of searching for the saddle point nearest to the local 
minimum of the potential energy E pot in normal coordinates (for more details, see [19, 
20]). It should be noted that the task of finding the saddle point in this case can be 
substantially simplified if we use the molecular dynamics data, because they provide 
information both about the isomer, into which one or another polycubane transforms 
upon the decomposition, and about the atomic configurations that arise in the course of 
this transformation. 



3. RESULTS AND DISCUSSION 

The binding energies of atoms in the polycubanes C4 + 4 n H 8 , which were calculated 
according to the formula [10] 

= (4 + An)E(C) + ggg) - £(C 4+4 nH 8 ) 
b 12 + An 1 ' 

are equal to 4.51, 4.59, 4.65, and 4.69 eV/atom for n — 2, 3, 4, and 5, respectively. For 
comparison, the binding energy in the CgH 8 cubane is Eb = 4.42 eV/atom [7]. Thus, 
the binding energy Eb increases with increasing n, which indicates an increase in the 
thermodynamic stability. As will be shown below, this is not accompanied by an increase 
in the kinetic stability. 

For each n = 2 - 5, we investigated the evolution of the C 4+4n H 8 cluster for 30 - 60 
different sets of initial velocities and displacements of atoms, which correspond to the 
temperature ranges T = 500 - 1700 K (n = 2), 300 - 1400 K (n = 3), 80 - 500 K (n = 4), 
and 50 - 200 K (n = 5). The temperature range every time was chosen in such a way that 
the lifetime r of the corresponding polycubane up to its decomposition would change in 



the interval from ~ 1 ps to ~ 0.1 /is, which corresponds to molecular dynamics steps in 
the range from ~ 3 x 10 3 to ~ 3 x 10 8 , respectively. The choice r ~ 1 ps as the lower 
boundary of this interval is motivated by the requirement r >> r , where r ~ 30 fs is 
the period of the fastest vibrations of the C4 + 4„H 8 clusters. The upper boundary of the 
aforementioned interval r ~ 0.1 fj,s was chosen rather arbitrarily and is dictated by the 
limited speed of the computers (for the C 8 H 8 cubane, it can be increased to ~ 1 /is [7]). 
As will be shown below, the significant difference in the temperature intervals (and in the 
mean values of temperatures in these intervals) at different values of n is associated with 
the different thermal stabilities of the C4 + 4„H 8 polycubanes studied in our work. 

Figure 4 presents the results of the direct numerical calculation of the lifetime r at 
different temperatures T for bicubane (n = 2), tricubane {n = 3), and tetracubane 
(n = 4). The data obtained for pentacubane (n = 5) do not present in this figure for 
the purpose of facilitating its perception. It can be seen that the dependences of lnr 
on 1/T for all values of n are well approximated by straight lines in accordance with the 
Arrhenius equation (1) without a correction for the finite sizes of the thermal reservoir [21- 
23]. The slopes of these lines determine the corresponding activation energies E a , whereas 
the points of their intersection with the ordinate axis determine the frequency factors A 
involved in formula (1). A statistical analysis of the data obtained from the "computer 
experiment" gives the following values of the aforementioned parameters: E a = 0.76 ±0.1, 
0.36 ± 0.05, 0.10 ± 0.02, and 0.07 ± 0.02 eV and A = io 14 - 79±0 - 22 , io 14 - 16±0 - 12 , Kp-osio.u 
i(]i3.72±o.72 s _1 for n — 2, 3, 4, and 5, respectively. Taking into account that E a = 
1.9 ±0.1 eV for the C 8 H 8 cubane [7], we can draw the conclusion that the activation energy 
E a of decomposition of the polycubanes in the C 4+4n H 8 family monotonically decreases 
with increasing n (Fig. 5). The observed decrease in the stability of the polycubanes is 



apparently associated with the fact that, beginning with n = 2 in the C 4+4n H 8 clusters, 
there arise purely carbon C-C-C bonds, which are bent at a right angle and do not contain 
(unlike the case of the C 8 H 8 cubane) hydrogen atoms. It can be seen from Fig. 4 that, for 
each value of n, the lifetime r exponentially decreases with an increase in the temperature 
T and that, at a fixed temperature T < 1000 K, it decreases (also exponentially) with an 
increase in the number n of C 8 cubes in the carbon skeleton of the polycubane. 

The frequency factor A with increasing number n also exhibits a tendency toward a 
decrease (according to the data reported in [7] for the C 8 H 8 cubane, the frequency factor 
is equal to io 16 03±0 - 36 s _1 ). Consequently, at a sufficiently high temperature T > 1000 
K, the lifetime of the polycubanes increases with an increase in their length, because, 
owing to the small ratio Ea/ksT, the preexponential factor in the Arrhenius equation (1) 
does not play a decisive role. Note, however, that interest expressed in the range of these 
temperatures is rather limited, because the absolute value of r in this case remains very 
small (~ 1 ps, Fig. 4). 

Let us analyze the structure and energy of the isomers formed upon the decomposition 
of the polycubanes. The molecular dynamics data indicate that the decomposition of 
the Ci2H 8 bicubane almost always occurs as a result of the synchronous breaking of two 
C-C bonds, which are directed perpendicular to the bicubane axis (and parallel to each 
other) and located at the center of the cluster. This results in the formation of an isomer 
shown in Fig. 6. Two opposite side faces of this isomer are hexagons formed by the 
C-C bonds, as is the case in nanotubes. The binding energy in the isomer increases to 
E b = 4.68 eV/atom; i.e., such isomer is energetically more favorable as compared to the 
bicubane (see formula (2)). In the further evolution, this isomer transforms into different 
quasi-two-dimensional clusters and no return to the initial configuration of the bicubane 



occurs. 

The decomposition of the Ci 6 H 8 tricubane begins with the sequential formation of the 
isomers depicted in Fig. 7. First, one pair of C-C bonds perpendicular to the tricubane 
axis (and parallel to each other) undergoes breaking, followed by the breaking of another 
pair of bonds that is perpendicular to the first pair and shifted relative to it along the 
tricubane axis by the length of the C-C bond. This leads to the formation of a cluster 
containing four hexagons composed of C-C bonds. The binding energy E b = 5.06 eV/atom 
of this cluster is higher than that of the tricubane. Then, this isomer transforms into 
different low-dimensional hydrocarbon complexes. We have never observed the recovery 
of the initial tricubane. 

The decomposition of the polycubanes C 4+4n H 8 with n > 3 occurs in a similar manner. 
Every time, after the synchronous breaking of pairs of transverse C-C bonds, the rectangles 
consisting of C-C bonds sequentially transform into hexagons so that, at n » 1, we 
eventually observe the formation of the (2,2) nanotube passivated with hydrogen atoms 
along the edges (see also [6]). 

Since the breaking of the first pair of C-C bonds leads to an irreversible decomposition 
of each polycubane, the height U of the energy barrier preventing this decomposition 
represents just the height of the barrier to the decomposition of this polycubane. Figure 
8 presents the results of the calculation of the potential energy E pot of the Ci 2 H 8 bicubane 
along the reaction coordinate connecting the initial metastable configuration with the 
configuration of the isomer shown in Fig. 6. The energy barrier height thus determined, 
i.e., U = 0.69 eV, coincides within the statistical error with the activation energy of 
thermally activated decomposition of the bicubane E a = 0.76 ±0.1 eV. Similarly, we have 
determined the height U = 0.39 eV of the energy barrier to the decomposition of the 



Ci 6 H 8 tricubane (Fig. 9). This value is also consistent with the corresponding activation 
energy E a = 0.36 ± 0.05 eV. It should be noted that, for the first isomer formed from 
the tricubane, there is a very shallow local minimum of the potential energy E pot with 
a depth of ~ 0.01 eV (Fig. 9). Consequently, in the molecular dynamics simulation of 
the evolution of the tricubane, we have observed configurations that are very close to this 
isomer for a very short time (~ 1 ps), followed by the formation of the second isomer 
Ci6H 8 (Fig. 7b). For tetracubane and pentacubane, we found the energy barrier heights 
U = 0.11 and 0.06 eV, respectively, which is also consistent with the activation energies 
E a = 0.10 ± 0.02 and 0.07 ± 0.02 eV, respectively. 

4. CONCLUSIONS 

Thus, the results presented in this paper on the numerical simulation of the dynamics of 
polycubanes C 4+ 4 n H 8 with different values of n indicate a rapid decrease in their thermal 
stability with an increase in the number n of cubic fragments C 8 in the carbon skeleton. 
For each member of the family of polycubanes with n = 2 - 5, we have determined the 
activation energy of its decomposition E a and the frequency factor A, which have allowed 
one to find the lifetimes r at an arbitrary temperature. In particular, at T = 300 K, 
we have r ~ 10 ms, 10 ns, 10 ps, and 0.1 ps for bicubane, tricubane, tetracubane, and 
pentacubane, respectively; i.e., an increase in the length of the polycubane leads to a rapid 
decrease in r. Therefore, it is hardly possible to conclude that, under normal conditions, 
there exist "nanotubes" with a square cross section (even if their edges are passivated 
with hydrogen). 

Meanwhile, for the lifetime of bicubane and tricubane at a temperature of 77 K, we 
have the estimates r ~ 10 35 and 10 9 s, respectively. This suggests that the possibility 



exists of monitoring these unique hydrocarbon clusters at cryogenic temperatures (here, 
we do not touch the problem of their synthesis, which is a separate and very difficult 
task). Of particular interest (by analogy with solid cubane s-C 8 H 8 ) is the preparation 
of anisotropic solid bicubane s-Ci2H 8 and solid tricubane s-Ci6H 8 , in which polycubanes 
oriented in one direction are bound to each other by means of van der Waals bonds. It 
is also interesting to consider the existence of a superbicubane and a supertricubane, in 
which structural elements are rectangular parallelepipeds C12 and Ci6, respectively. 
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Fig. 1. CgHg cubane 



(1) carbon atoms and (2) hydrogen atoms. 




Fig. 2. A fragment of the supercubane. 




Fig. 3. (a) Bicubane C^Hg and (b) tricubane CieHg. The notation of the atoms is the 
same as that used in Fig. 1. 
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Fig. 4. Dependences of the logarithm of the lifetime r for the (1) bicubane Ci 2 H 8 , (2) 
tricubane Ci 6 H 8 , and (3) tetracubane C 2 oH 8 on the reciprocal of the initial temperature 
T. Solid lines are the corresponding linear approximations obtained by the least-squares 
method. 




Fig. 5. Dependence of the activation energy of decomposition of the C4 + 4„H 8 cluster 
on the number n of C 8 cubes in the carbon skeleton (data for n — 1 are taken from [7]). 




Fig. 6. Isomer Ci^Hg into which the bicubane transforms upon the decomposition. The 
notation of the atoms is the same as that used in Fig. 1. 



Fig. 7. Isomers CieHg into which the tricubane sequentially transforms upon the 
decomposition: (a) metastable isomer and (b) low-energy isomer. The notation of the 
atoms is the same as that used in Fig. 1. 
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Fig. 8. Potential energy E pot of the Ci2H 8 cluster in the vicinity of the configuration 
of the bicubane: (1) bicubane (Fig. 3a), (2) saddle point that determines the barrier 
U = 0.69 eV for the decomposition of the bicubane, and (3) low-energy isomer (Fig. 6). 
The reaction coordinate is chosen as the half the sum of increments in the lengths of two 
C-C bonds broken during the decomposition of the bicubane. 
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Fig. 9. Potential energy E pot of the Ci 6 H 8 cluster in the vicinity of the configuration 
of the tricubane: (1) tricubane (Fig. 3b), (2) saddle point that determines the barrier 
U = 0.39 eV for the decomposition of the tricubane, (3) metastable isomer (Fig. 7a), 
(4) saddle point that determines the barrier U = 0.008 eV for the decomposition of this 
isomer, and (5) low-energy isomer (Fig. 7b). The reaction coordinate is chosen as the 
half the sum of increments in the lengths of four C-C bonds sequentially broken during 
the decomposition of the tricubane (see text). 



